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1.  Introduction 


The  archiving  of  ionospheric  characteristics  at  uneven  data 
rates  was  presented  by  Gamache  and  Reinisch  [1989a]  at  the 
International  Workshop  on  Digital  lonogram  Data  Formats  for  World 
Data  Center  Archiving  [Gamache  and  Reinisch,  1989b]  as  part  of  the 
URSI  Working  Group  G.4:  Ionospheric  Informatics.  Because  of  the 
uneven  time  sampling,  the  database  has  records  of  varying  sizes. 
Historically,  this  has  created  problems  in  both  the  creation  of  and  in 
later  utilization  of  the  files.  The  IIWG  proposed  database  avoids 
these  problems  by  having  the  header  records  serve  as  a  key  to 
encoding  the  remainder  of  the  data.  The  database  is  standard  ASCII 
characters  and  is  created,  serviced,  and  used  via  FORTRAN  programs. 
For  a  full  description  of  the  database  and  the  philosophy  behind  its 
creation  see  Gamache  and  Reinisch  1898a.  Table  1  gives  the  revised 
IIWG  proposed  structure  of  the  database. 

2.  Description  of  characteristics  database 

Before  describing  the  structure  of  the  database,  it  is  important 
to  give  some  definitions  to  aid  in  the  discussion.  A  database  record, 
DR,  is  defined  as  the  collection  of  all  characteristics  data  for  a  month. 
The  database  is  created  by  a  FORTRAN  program  (this  can  also  be 
done  in  C  or  any  high  level  language)  and  each  FORTRAN  WRITE 
statement  creates  a  FORTRAN  record,  FR,  from  which  the  database  is 
comprised.  For  ease  of  use,  transportability,  editing,  and  printing, 
the  FR  length  was  set  to  120  characters  maximum.  When  a  DR  is 
created  on  a  hard,  floppy,  or  optical  disk  it  is  called  a  file  and  is  a 
single  entity.  When  created  on  a  magnetic  tape,  the  file  is 
constructed  in  a  number  of  blocks.  Tape  handling  routines  put  a  0.6 
inch  gap  between  blocks.  For  magnetic  tape  the  block  size  has  been 
set  by  adjusting  the  block  size  to  maximize  the  ratio  of  the  inter¬ 
block  gaps  and  data  stored.  On  some  systems  the  block  size  is 
governed  by  software  that  limits  the  counter  to  a  four  byte  word 
yielding  a  maximum  block  size  of  9999  bytes.  For  transportability,  a 
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Table  1.  URSI  IIWG  Database  Structure  for  Flexible  Data  Rates 


Data 

Gioup 

FORTRAN 

Record# 

Format 

Description 

1 

A30 

Station  Name 

1 

A5 

Station  code 

1 

14 

Meridian  time  used  by  station  on  records 

! 

1 

F5.1 

Latitude  N 

1 

F5.1 

Longitude  E 

1 

AlO 

Scaling  type:  Manual /Automatic 

1 

AlO 

Data  editing:  Edited/Non-edited/Mixed 

1 

A30 

lonosonde  system  name 

2,3 

3014* 

Year 

2 

4,i 

12A10* 

Month 

Number  of  days  in  the  month,  M 

Number  of  Characteristics 

Numbers  of  measurements  total 

Numbers  of  measurements  for  each  of 
the  M  days,  Nm 

List  of  characteristics 

i+l,j 

12A10* 

Dimensions 

j+l,k 

60A2* 

List  of  corresponding  URSI  codes 

k+1,1 

20(312)* 

The  Nm  sample  times  Hh:Mm:Ss 

3 

• 

•  • 

for  each  of  the  M  days 

l+l,m 

24(13,2  AD* 

The  Ni  values  of  characteristic  1  for  day  1 

• 

•  • 

•••repeated  for  each  of  the  M  days 

m+1 

24(I3,2AD 

Hourly  Medians  for  characteristic  1 

m+2 

24(12,2  AD 

The  counts  for  the  hourly  medians,Range 

m+3 

24(13,2  AD 

Upper  quartile 

m44 

24(I3,2AD 

Lower  quartile 

m+5 

24(I3,2AD 

Upper  decile 

m+6 

24(13,2AD 

Lower  decile 

m+7,n 

24(I3,2AD* 

The  N2  values  of  characteristic  2  for  day  1 

• 

• 

•  •  etc. 

• 

•  • 

repeated  for  each  characteristic 

•  •  • 

*  Format  is  repeated  as  many  times  as  necessary  to  read/write  the  data. 
Revised:!!  April  !99!. 
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block  size  of  9000  bytes  was  chosen.  This,  for  example,  is  capable  of 
storing  two  characteristics  of  routine  hourly  measurements  within  a 
block  on  magnetic  tape. 

The  structure  of  the  IIWG  characteristics  database  for  a 
database  record  is  displayed  in  Table  1.  A  given  DR  is  comprised  of  a 
number  of  FRs,  the  total  number  of  FRs  being  determined  by  the 
number  of  days  in  the  month,  the  number  of  measurements  made 
each  day,  and  the  number  of  characteristics  being  archived.  The 
database  record  is  constructed  with  the  data  being  in  three  groups. 
The  first  group  is  informative  and  encoding  data.  This  starts  with  a 
FR  containing  the  Station  Name,  (A30)  format,  where  the  data  were 
recorded  and  the  Station  Code  in  (AS)  format;  the  meridian  time 
used  by  the  station  to  indicate  if  time  is  recorded  in  UT  or  LT  on  the 
records  is  given  next  in  (14)  format  followed  by  the  station 
coordinates.  Latitude  N  and  Longitude  E  both  in  (F5.1)  format; 
next  there  are  two  (AlO)  format  variables  describing  the  Scaling 
type,  this  takes  the  value  Manual  or  Automatic,  and  the  Data 
editing  variable  which  can  be  Edited,  Non-Edited,  or  Mixed.-,  last  in 
this  FR  is  space  for  the  lonosonde  system  name  in  (A30)  format. 

FORTRAN  records  two  and  three  contain  most  of  the 
information  for  encoding  the  DR.  A  repeating  (3014)  format  is  used 
to  write  both  records.  They  start  with  the  Year,  Month,  Number  of 
days  in  the  month  (NN),  the  Number  of  characteristics  (Nchar), 
archived  in  this  DR,  the  Total  number  of  measurements  reported 
for  the  DR,  and  the  Number  of  daily  measurements  made  for 
each  of  the  NN  days.  FR  number  four  to  record  i  use  a  repeating 
format  of  (12 AlO)  to  list  the  Names  of  the  particular  characteristics 
being  archived.  There  are  Nchar  of  them,  hence  the  FRs  are  repeated 
as  needed.  For  example  if  one  were  archiving  only  the  critical 
frequencies  foF2,  foFl,  and  foE,  Nchar  would  be  three,  and  the 
characteristics  list  would  be  ’  foF2'  ’  foFl'  '  foE'.  A  list  of 
the  names  of  the  characteristics,  the  units,  and  URSI  codes  taken 
from  UAG23  (Piggott  and  Rawer,  1978)  are  given  in  Table  2;  The 
URSI  list  has  been  enhanced  with  characteristics  that  are  scaled  by 
the  Digisonde  ARTIST  [Tang  et  al.  1989].  FORTRAN  records  i  +  1  \o  k 
give  the  Dimensions  corresponding  to  the  characteristics  list  (see 
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Table  2),  these  are  in  (12A10)  repeating  format.  The  last  FRs  in  the 
first  group  are  for  the  URSI  codes  specified  for  each  of  the 
characteristics  (see  Table  2)  and  are  written  in  (60A2)  repeating 
format. 

From  the  information  in  this  first  group  one  knows 
immediately  how  many  data  for  the  time  or  for  each  characteristic 
are  to  be  read.  From  the  number  of  measurements  for  each  day  the 
time  data  can  be  separated  into  the  times  for  the  individual  days  of 
the  month  and  the  measured  characteristics  can  uniquely  be 
associated  with  a  given  time  on  a  given  day. 

The  second  group  of  FRs  contain  the  measurement  times  for  the 
month.  With  uneven  time  spacing  the  measurement  times  must  be 
recorded  to  associate  with  the  reported  characteristics.  This  requires 
that  hours,  minutes,  and  seconds  of  each  measurement  be  entered 
into  the  database.  To  conserve  space,  the  times  are  written  once  per 
month  and  the  reported  characteristics  are  written  to  correspond  to 
these  times.  FORTRAN  records  k+1  to  /  are  used  for  the  measurement 
times.  They  are  written  in  a  (30(312))  repeating  format 
corresponding  to  the  hours,  minutes,  and  seconds,  HHMMSS,  of  the 
measurements.  The  number  of  FRs  needed  for  this  is  determined  by 
the  data  sampling  rate  for  the  month. 

The  third  group  of  FRs  contains  the  actual  values  of  the 
characteristics  and  the  corresponding  hourly  medians  and  statistics. 
The  group  is  comprised  of  a  number  of  FRs  for  each  archived 
characteristic  which  are  repeated  for  each  characteristic.  The  order 
of  the  characteristics  follows  that  given  in  the  "List  of 
characteristics".  On  a  per  characteristic  basis,  for  each  characteristic 
one  has  the  Ni  values  of  the  characteristic  for  day  1  corresponding  to 
the  reported  measurement  times  for  day  1.  These  are  followed  by 
the  values  for  day  2,  day  3,  ...  for  each  of  the  NN  days  of  the  month. 
The  characteristics  are  written  in  a  repeating  (24(13 ,2 A1 ))  format 
corresponding  to  the  integer  value  (13)  of  the  characteristic  and  the 
qualifying  and  descriptive  letters  [see  UAG  23]. 

The  IIWG  Workshop  suggested  the  use  of  two  slashes,  //,  in 
place  of  the  qualifying  and  descriptive  letters  for  monthly 
characteristics  data  that  were  autoscaled  but  not  checked  or  "edited". 
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Table  2.  List  of  Characteristics,  URSI  codes,  and  Dimensions 


GROUP  CHARACTERISTIC  DIMENSIW  REFERENCE  DEFINITION 

ARTIST  URSI  UAG23 

Name  #  Name  # 


FI 


1 

foF2 

00 

.1 

MHz 

1.11 

The  ordinary  wave  critical  frequency  of  the 
highest  stratification  in  the  F  region 

fxF2 

01 

.1 

MHz 

1.1 1 

The  extraordinary  wave  critical  frequency 

fzF2 

02 

.1 

MHz 

1.11 

The  z-mode  wave  critical  frequency 

M(D) 

3 

M3000F2 

03 

.01 

MHz 

1.50 

The  maximum  usable  frequency  at  a  defined 
distance  divided  by  the  critical  frequency  of 
that  layer 

hpF2 

12 

h'F2 

04 

km 

1.33 

The  minimum  virtual  height  of  the  ordinary 
wave  trace  for  the  highest  stable  stratification 
in  the  F  region 

hpF2 

05 

km 

1.41 

The  virtual  height  of  the  ordinary  wave  mode 
at  the  frequency  given  by  0.834  of  foF2  (or 
other  7.34) 

h’Ox 

06 

km 

1.39 

The  virtual  height  of  the  x  trace  at  foF2 

MUF(D) 

4 

MUF3000F2 

07 

.1 

MHz 

1.5C 

The  standard  transmission  curve  for  3000  km 

he 

08 

km 

1.42 

The  height  of  the  maximum  obtained  by 
fitting  a  theoretical  h'F  curve  for  the  parabola 
of  best  fit  to  the  observed  ordinary  wave  trace 
near  foF2  and  correcting  for  underlying 
ionization 

qc 

09 

km 

7.34 

Scale  height 

2 

foFl 

10 

.01 

MHz 

1.13 

The  ordinary  wave  FI  critical  frequency 

fxFl 

11 

.01 

MHz 

1.13 

The  extraordinary  wave  FI  critical  frequency 

12 

not  used 

M3000F1 

13 

.01 

MHz 

1.50 

See  Code  03 

h'Fl 

14 

km 

1.30 

The  minimum  virtual  height  of  reflection  at  a 
point  where  the  trace  is  horizontal 

15 

not  used 

hpF 

11 

h'F 

16 

km 

1.32 

The  minimum  virtual  height  of  the  ordinary 
wave  trace  taken  as  a  whole 

MUF3000F1 

17 

.1 

MHz 

1.5C 

See  Code  07 

18 

not  used 

19 

not  used 

9 

foE 

20 

.01 

MHz 

1.14 

The  ordinary  wave  critical  frequency  of  the 
lowest  thick  layer  which  causes  a 
discontinuity 

21 

not  used 

foE2 

22 

.01 

MHz 

1.16 

The  critical  frequency  of  an  occulting  thick 
layer  which  sometimes  appears  between  the 
normal  E  and  FI  layers 

23 

not  used 

hpE  13 

h'E 

24 

km 

1.34 

The  minimum  virtual  height  of  the  normal  E 
layer 

25 

not  used 

h'E2 

26 

km 

1.36 

The  minimum  virtual  height  of  the  E2  layer 

trace 

27 

not  used 

28 

not  used 

29 

not  used 
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Table  2. .continued 


GROUP  CHARACTERISTIC  DIMENSKM4 


ARTIST 

URSI 

Name  # 

Name 

# 

Es 

6 

foEs 

30 

.1  MHz 

fxEs 

31 

.1  MHz 

fbEs 

32 

.1  MHz 

ftEs 

33 

.1  MHz 

hpEs  1 4 

h'Es 

Type  Es 

34 

35 

36 

37 

38 

39 

km 

Other 

foFl.5 

40 

41 

.01  MHz 

5 

fmin 

42 

.1  MHz 

M3000F1.5 

43 

.01  MHz 

h'Fl.5 

44 

45 

46 

km 

fm2 

47 

.1  MHz 

hm 

48 

km 

fm3 

49 

.1  MHz 

Spread  F/Oblique 

fol 

50 

.1  MHz 

10 

fxl 

51 

.1  MHz 

frnl 

52 

.1  MHz 

M30001 

53 

.01  MHz 

hi 

54 

55 

56 

km 

dfs 

57 

58 

59 

.1  MHz 

REFERENCE  DEFINITION 
UAG23 


1.17  The  highest  ordinary  wave  frequency  at  which 
a  mainly  continuous  Es  trace  is  observed 

1.17  The  highest  extraordinary  wave  frequency  at 
which  a  mainly  continuous  Es  trace  is 
observed 

1.18  The  blanketing  frequency  of  the  Es  layer 
Top  frequency  Es  any  mode. 

1.35  The  minimum  height  of  the  trace  used  to  give 
foEs 
not  used 

7.26  A  characterization  of  the  shape  of  the  Es  trace 
not  used 
not  used 
not  used 


1.12  The  ordinary  wave  critical  frequency  of  the 
intermediate  stratification  between  FI  and  F2 
not  used 

1.19  The  lowest  frequency  at  which  echo  traces  are 
observed  on  the  ionogram 

1.50  See  Code  03 

1.38  The  minimum  virtual  height  of  the  ordinary 
wave  trace  between  foFl  and  foFl.5  (equals 
h'F2  7.34) 
not  used 
not  used 

1.14  The  minimum  frequency  of  the  second  order 
trace 

7.34  The  height  of  the  maximum  density  of  the  F2 
layer  calculated  by  the  Titheridge  method 

1 .25  The  minimum  frequency  of  the  third  order 
trace 


1 .26  The  lop  ordinary  wave  frequency  of  spread  F 
traces 

1.21  The  lop  frequency  of  spread  F  traces 

1.23  The  lowest  frequency  of  spread  F  traces 
1 .50  See  Code  03 

1.37  The  minimum  slant  range  of  the  spread  F 
traces 
not  used 
not  used 

1.22  The  frequency  spread  of  the  scatter  pattern 

7.34  Frequency  range  of  spread  fxI-foF2 

not  used 
not  used 
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Table  2..continued 


GROUP  CHARACTERISTIC 

DIMENSIW 

REFERENCE  DEFINITION 

ARTIST 

URSI 

UAG23 

Name  # 

Name 

# 

N(h) 

30 

fh’F2 

60 

.1  MHz 

7.34 

The  frequency  at  which  h'F2  is  measured 

29 

fh’F 

61 

.1  MHz 

7.34 

The  frequency  at  which  h'F  is  measured 

62 

not  used 

h'mFl 

63 

km 

7.34 

The  maximum  virtual  height  in  the  o-mode  FI 
cusp 

hi 

64 

km 

7.34 

True  height  at  fl  Titheridge  method 

h2 

65 

km 

7.34 

True  height  at  f2  Titheridge  method 

h3 

66 

km 

7.34 

True  height  at  D  Titheridge  method 

h4 

67 

km 

7.34 

Tiue  height  at  f4  Titheridge  method 

h5 

68 

km 

7.34 

True  height  at  f5  Titheri'^  e  method 

H 

69 

km 

7.34 

Effective  scale  height  at  hmF2  Titheridge 
method 

T.E.C. 

12000 

70 

e/cm^ 

7.34 

Ionospheric  electron  content  Faraday 
technique 

I 

71 

e/cm^ 

7.34 

Total  electron  content  to  geostationary 
satellite 

Ixxxx 

72 

e/cm^ 

7.34 

Ionospheric  electron  content  to  height  xxxx 

73 

not  used 

74 

not  used 

75 

not  used 

76 

not  used 

77 

not  used 

78 

not  used 

T 

79 

e/cm^ 

7.34 

Total  sub-peak  content  Titheridge  method 

AUTOMATIC 

7 

FMINF 

80 

.01  MHz 

Minimum  frequency  of  F  trace  (50  kHz 
increments)  Equals  fbEs  when  E  present 

8 

FMINE 

81 

.01  MHz 

Minimum  frequency  of  E  trace  (50  kHz 
increments). 

15 

HOM 

82 

km 

Parabolic  E  region  peak  height 

16 

YM 

83 

km 

Parabolic  E  region  semi-thickness 

17 

QF 

84 

km 

Average  range  spread  of  F  trace 

18 

QE 

85 

km 

Average  range  spread  of  E  trace 

22 

FF 

86 

.01  MHz 

Frequency  spread  between  fxF2  and  f.xl 

23 

FE 

87 

.01  MHz 

As  FF  but  considered  beyond  foE 

25 

fMUF3000 

88 

.01  MHz 

MUF(D)/obliquity  factor 

26 

h'MUF3000 

89 

km 

Virtual  height  at  fMUF 
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i.e.  where  no  quality  control  procedure  has  been  applied.  This  code 
has  been  extended  to  consider  data  that  have  been  edited  but  no 

descriptive  or  qualifying  letters  introduced.  With  two  positions  to  fill 

and  the  use  of  a  single  or  double  slash  there  are  four  codes  which  can 
be  defined.  The  first  is  no  slashes  implying  the  use  of  the  descriptive 
or  qualifying  letters.  The  next  is  the  use  of  two  slashes  which 
signifies  no  editing.  The  third  choice  is  to  put  a  slash  in  the  first 
position  followed  by  a  blank.  This  is  used  to  signify  autoscaled  data 
that  have  been  edited  but  no  descriptive  or  qualifying  letters  are 
used.  The  last  possibility  is  a  blank  in  the  first  position  followed  by 

the  slash.  This  is  not  currently  used  thus  it  leaves  the  possibility  for 

future  extension  of  the  code  The  codes  are  summarized  in  Table  3. 

Table  3.  IIWG  Codes  for  the  Descriptive  and  Qualifying  Fields  of 
the  Characteristics. 


Symbolic  Code _ Description 


Q  D 

Qualifying  and  descriptive  letters  used  according  to 
UAC,  #23 A. 

1  _ 

Autoscaled  data,  edited  but  no  qualifying  and 
descriptive  letters  used. 

-1 

No  current  meaning,  for  future  extension. 

1  d 

Autoscaled  data,  no  editing,  no  qualifying  and 
descriotive  letters  used. 

The  actual  values  of  the  characteristics  can  be  obtained  by 
multiplying  the  integer  value  by  the  value  found  in  the 
corresponding  "dimension  list”  (group  2)  of  the  database  record  (see 
Table  2).  Thus  a  value  of  86  reported  for  foF2  is  multiplied  by  the 
dimension  0.1  MHz  to  give  a  foF2  value  of  8.6  MHz. 

Immediately  following  the  characteristics  data  are  the  hourly 
medians  given  in  a  (24(13 ,2 A1 ))  format;  the  counts  for  the  hourly 
medians  and  the  range  in  (24(12 ,13))  format;  the  upper  quartiles 
in  a  (24(13 ,2  A1 ))  format;  the  lower  quartiles  in  a  (24(13 ,2A1 )) 
format;  then  the  upper  deciles  in 'd  (24(13 ,2AI ))  format;  and  finally 
the  lower  quartiles  again  in  a  (24(13 ,2A1 ))  format. 
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The  above  FRs  are  repeated  for  each  characteristic  given  in  the 
"characteristics  list."  This  completes  the  DR,  i.e.  a  month  of 
characteristics  data. 

The  reading  and  writing  of  the  database  records  to  magnetic 
tape  or  other  mediums  are  accomplished  by  FORTRAN  (or  other  high 
level  languages)  programs  distributed  with  the  database  (see  the 
Appendices).  The  uneven  time  spacing  leads  to  DR  lengths  that 
change  with  the  data,  but  as  shown  this  is  no  problem  given  the 
structure  of  the  database.  This  is  an  important  factor  when 
transportability  is  considered. 

The  advantages  of  this  structure  are:  the  flexibility  of  the 
resulting  database  allows  data  rates  to  change  from  month  to  month 
and/or  from  station  to  station.  The  database  is  ANSII  standard  and 
transportable  and  is  easily  maintained  and  utilized  by  FORTRAN 
routines  provided.  The  structure  is  very  flcAible  and  most 

importantly  does  not  lead  to  data  being  discarded. 

3.  Implementation  for  ADEP  data 

The  creation  of  the  characteristics  database  has  been 
automated  for  ULCAR  ADEP  data.  The  ADEP  data,  a  structured 
database  developed  by  Gamache  et.  al.  [1989c],  was  also  presented  at 
the  Lowell  Workshop.  At  the  workshop  it  was  suggested  that  this 
database  be  adopted  as  a  standard  for  archiving  both  manual  and 
autoscaled  ionogram  data.  To  accomplish  this  several  enhancements 
were  added  to  the  structure.  In  Table  4,  the  suggested  structure  for 
the  IIWG  ADEP  database  is  given;  this  structure  has  been 

programmed  into  the  ADEP  program.  This  is  the  third  change  in  the 
structure  of  the  ADEP  database,  this  however  has  not  been  a  problem 
because  the  database  was  created  to  be  upwards  compatible.  Several 
programs  exist  to  convert  the  original  ADEP  data  (type  1,  1988)  to 
the  Standard  ADEP  Output  (SAO)  data  (type  2,  1989),  and  now  from 
the  SAO  format  to  the  IIWG  recommended  format  (type  3,  1990), 
without  loss  of  any  data.  The  Workshop’s  recommendation  was  that 
the  ADEP  structure  for  scaled  ionogram  data  serve  as  a  guide  for 
other  ionosonde  systems. 
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Table  4.  IIWG  ADEP  Structure  as  of  January  31,  1990 


Code 

Format 

Description 

2(4013) 

DATA  HLE  INDEX 

1 

N(16F7.3) 

GEOPHYSICAL  CONSTANTS 

2 

N(A120) 

SYSTEM  DESCRIPTION:  Version  numbers,  etc. 

3 

N(120Z1) 

lONOGRAM  SOUNDING  SETTINGS  (PREFACE) 

4 

N(15F8.3) 

SCALED  IONOSPHERIC  PARAMETERS 

5 

N(60I2) 

ANALYSIS  FLAGS 

6 

N(16F7.3) 

DOPPLER  TRANSLATION  TABLE 

0-TRACE  POINTS  -  F2  LAYER 

7 

N(40I3) 

VIRTUAL  HEIGHTS 

8 

N(40I3) 

TRUE  HEIGHTS 

9 

N(60I2) 

AMPLITUDES 

10 

N(120I1) 

DOPPLER  NUMBER 

11 

N(20F6.3) 

FREQUENCY  TABLE 

O-TRACE  POINTS  -  FI  LAYER 

12 

N(40I3) 

VIRTUAL  HEIGHTS 

13 

N(40I3) 

TRUE  HEIGHTS 

14 

N(60I2) 

AMPLITUDES 

15 

N(120I1) 

DOPPLER  NUMBER 

16 

N(20F6.3) 

FREQUENCY  TABLE 

O-TRACE  POINTS  -  E  LAYER 

17 

N(40I3) 

VIRTUAL  HEIGHTS 

18 

N(40I3) 

TRUE  HEIGHTS 

19 

N(60I2) 

AMPLITUDES 

20 

N(120I1) 

DOPPLER  NUMBER 

21 

N(20F6.3) 

FREQUENCY  TABLE 

X-TRACE  POINTS  -  F2  LAYER 

22 

N(40I3) 

VIRTUAL  HEIGHTS 

23 

N(60f2) 

AMPLITUDES 

24 

N(120I1) 

DOPPLER  NUMBER 

25 

N{20F6.3) 

FREQUENCY  TABLE 

X-TRACE  POINTS  -  FI  LAYER 

26 

N(40I3) 

VIRTUAL  HEIGHTS 

27 

N(60I2) 

AMPLITUDES 

28 

N(120n) 

DOPPLER  NUMBER 

29 

N(20F6.3) 

FREQUENCY  TABLE 

X-TRACE  POINTS  -  E  LAYER 

30 

N(40I3) 

VIRTUAL  HEIGHTS 

31 

N(60I2) 

AMPLITUDES 

32 

N(120I1) 

DOPPLER  NUMBER 

33 

N(20F6.3) 

FREQUENCY  TABLE 

34 

N(60I2) 

MEDIAN  AMPLITUDE  OF  F  ECHO 

35 

N(6012) 

MEDIAN  AMPLITUDE  OF  E  ECHO 

36 

N(60I2) 

MEDIAN  AMPLITUDE  OF  ES  ECHO 

37 

N(13E9.4E1) 

TRUE  HEIGHT  F2  LAYER  COEFFICIENTS 

38 

N(13E9.4E1) 

TRUE  HEIGHT  FI  LAYER  COEFFICIENTS 

39 

N(13E9.4E1) 

TRUE  HEIGHT  E  LAYER  COEFFICIENTS 

40 

N(13E9.4E1) 

VALLEY  DESCRIPTION 

41 

N(120I1) 

EDIT  FLAGS 
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The  monthly  characteristics  program  was  developed  to  read  in 
one  month  of  ADEP  data  and  create  the  IIWG  characteristics 
database.  This  is  accomplished  by  reading  data  for  one  ionogram  at  a 
time  and  storing  the  measurement  time  and  characteristics  data  in 
arrays.  The  program  keeps  track  of  the  number  of  measurements 
made  each  day  of  the  month,  which  changes  with  the  measurement 
schedule.  Currently  the  arrays  in  the  program  for  this  are 
dimensioned  31X300  supporting  a  full  month  at  a  5  minute  schedule 
(288  measurements/day).  If  more  measurements  are  made  the 
array  boundaries  must  be  increased  accordingly. 

The  locations  of  the  hourly  measurements  (hours  0  to  23)  in 
the  measurements/day  list  are  found  and  stored  for  use  in  the 
median  routine.  This  array  keys  the  locations  of  the  hourly 
characteristics  in  the  larger  monthly  characteristics  arrays. 
Programming  in  this  fashion  is  very  memory  intensive,  however  it  is 
logical  and  produces  a  robust  and  error  free  algorithm.  For  example, 
the  data  supplied  to  the  algorithm  need  not  be  ordered  in  time  for 
the  algorithm  to  work  properly.  This  is  important  because  a  month 
of  ADEP  data,  several  megabytes,  seldom  exist  as  a  file  complete  for 
the  month  but  must  be  collated  into  a  single  file.  Errors  made  in 
collating  do  not  affect  the  characteristics  program.  Thus  the  results 
are  independent  of  how  the  data  for  the  month  were  put  together. 

The  number  of  FRs  necessary  for  archiving  each  characteristic 
is  variable  from  month  to  month  depending  on  the  measurement 
schedule.  The  program  makes  all  the  necessary  provisions  and  the 
key  information  is  written  into  the  first  group  of  FRs  for  future 
encoding.  In  this  way  relatively  little  space  is  used  for  null  data 
compared  to  the  previous  URSI  characteristics  files  [Rodger  1984], 

The  program  CHARS.FOR  was  written  in  FORTRAN??  and  will 
run  on  any  system  supporting  this  language.  It  was  developed  on  an 
IBM-PC/AT  compatible  system.  Because  of  the  large  arrays  in  the 
program,  COMMON  blocks  were  constructed  in  groups  each  with  less 
than  64K  requirements.  Data  not  in  common  blocks  do  not  exceed 
64K,  however  for  saftey  the  program  was  compiled  using  Microsoft 
FORTRAN  option  /Gt3000  which  partitions  this  data  so  that  the  64K 


constraint  is  never  exceeded.  It  is  recommended  that  similar  options 
be  used  on  other  systems  when  necessary. 

The  program  was  tested  using  several  months  of  simulated 
ADEP  data.  The  simulated  data  was  constructed  so  that  the  value 
assigned  to  a  characteristic  was  related  to  the  measurement  time. 
For  example,  the  critical  frequency  of  the  F2  layer  is  defined  by 

foF2  =  Hour  +  2.  +  Day/10.  (I) 

where  Hour  is  the  hour  (0  to  23)  and  Day  is  the  day  of  the  month 
(1  to  31).  This  allows  the  resulting  characteristics  files  and  median, 
quartile,  etc.  results  to  be  thoroughly  understood  and  checked.  The 
program  was  tested  with  input  files  that  varied  the  number  of 
measurements  per  day  from  24  to  288.  In  the  tests,  the 
characteristics  and  the  medians,  quartiles,  deciles,  counts  and  ranges 
were  checked  and  found  to  be  correct. 

In  additional  tests,  the  days  of  the  month  were  scrambled  and 
the  characteristics  files  generated  again.  The  order  of  the 

measurement  times  and  characteristics  were  unchanged  compared  to 
the  previous  file  and  the  medians,  quartiles,  etc.  were  unchanged  by 
the  reordering  of  the  days.  A  similar  test  was  done  where  in 
addition  to  scrambling  days,  the  measurement  times  within  each  day 
were  scrambled.  The  measurement  times  and  characteristics  in  the 
measurements/day  lists  had  a  different  order  but  they  had  a  one  to 
one  correspondences  in  the  FRs.  Of  course  the  hourly  measurement 
indication  vectors  were  very  different  but  the  medians,  quartiles, 
counts,  etc.  were  unchanged  for  the  month.  It  is  not  suggested  that 
ADEP  data  files  be  created  without  time  ordering,  however  if  this 
does  happen  it  will  not  affect  the  final  database  record. 

Appropriate  coding  has  been  added  as  edit  flag  information  to 
the  IIWG  ADEP  format.  Table  4  #41.  For  these  flags  a  0  implies 
autoscaling  and  a  1  indicates  the  corresponding  data  has  been 
manually  edited.  The  list  of  data  for  the  edit  flags  is  shown  in  Table 
5.  The  flag  position  is  from  the  position  of  the  scaled  characteristic  in 
the  ARTIST  and  ADEP  list  with  the  addition  of  flags  for  the  traces  and 
true  height.  A  program  was  written  to  take  the  standard  ADEP 
output  files,  type  2,  and  add  this  information  to  the  database  to 


TABLE  5.  ADEP  edit  flags. 


Edit  flag 
Position 

Description 

1 

foF2 

F2  layer  critical  frequency  calculated  by 
hyperbolic  fit 

2 

foFl 

FI  layer  critical  frequency 

3 

M(D) 

M-factor,  MUF(D)/foF2,  for  distance  D 

4 

MUF(D) 

Maximum  usable  frequency  for  distance  D 

5 

fmin 

Minimum  frequency  for  E  or  F  echoes 

6 

foEs 

Es  layer  critical  frequency 

7 

fminF 

Minimum  frequency  of  F-trace 

8 

fminE 

Minimum  frequency  of  E-trace 

9 

foE 

E  layer  critical  frequency 

10 

fxl 

Maximum  frequency  of  F-trace 

11 

h’F 

Minimum  virtual  height  of  F  trace 

12 

h’F2 

Minimum  virtual  height  of  F2  trace 

13 

h’E 

Minimum  virtual  height  of  E  trace 

14 

h'Es 

Minimum  virtual  height  of  Es  layer 

15 

HOM 

Peak  of  E  layer  using  parabolic  model 

16 

Ym 

Corresponding  half  thickness  of  E  layer 

17 

CF 

Average  range  spread  of  F-trace 

18 

CE 

Average  range  spread  of  E-trace 

19 

Down  F2 

Lowering  of  F-trace  maximum  to  leading 
edge 

20 

Down  E 

Lowering  of  E-trace  maximum  to  leading 
edge 

21 

Down  Es 

Lowering  of  Es-trace  maximum  to  leading 
edge 

22 

FF 

Frequency  spread  between  fxF2  and  fxl 

23 

FE 

As  FF  but  considered  beyond  foE 

24 

D 

Distance  used  for  MUF  calculation t 

25 

fMUF(D) 

MUF(D)/obliquity  factor* 

26 

h'MUF(D) 

Virtual  height  at  fMUF 

27 

foF2c 

correction  to  add  to  foF2  to  get  actual  foF2 

28 

foEp 

Predicted  foE 

29 

f(h’F) 

Frequency  at  which  hminF  occurs 

30 

f(h'F2) 

Frequency  at  which  hminF2  occurs 

31 

foFlp 

Predicted  foFl 

32 

Zmax 

Peak  height  F2  layer 

50 

F2  trace 

F2  trace  points  were  edited 

51 

FI  trace 

FI  trace  points  were  edited 

52 

E  trace 

E  trace  points  were  edited 

53 

true  height 

true  height  was  recalculated  with  edited  traces 

54 

Es  trace 

Es  trace  points  were  edited 

t  Normally  3000  km 

*  Obliquity  factor(h',  D)  is  the  ratio  of  frequencies  for  vertical  and  oblique 
propagation  to  distance  D  with  virtual  height  h’. 


produce  the  type  3  files,  IIWG  ADEP  database  files.  With  this  as 
input  the  addition  of  slashes  when  necessary  is  fully  automated. 

We  have  used  the  term  "autoscaled  data"  in  the  discussion 
since  our  experience  is  based  on  the  Digisonde  ARTIST-  scaled 
ionogram  data.  All  definitions  apply,  however,  as  well  to  manually 
or  semi-automatically  digitized  data. 

4.  Supplemental  Programs 

In  addition  to  the  FORTRAN  program  that  creates  the  database 
records,  an  additional  program  was  written  to  read  the  database 
records  and  write  out  the  data  along  with  some  descriptive 
information.  This  program  can  also  be  used  as  a  seed  program  to 
encode  the  characteristics  database  for  other  applications  and  studies 
involving  the  database  files,  allowing  lesser  experienced 
programmers  quick  access  to  the  data.  Like  the  database,  the  output 
files  are  written  with  a  maximum  record  length  of  120  characters 
allowing  a  study  of  the  data.  Currently  fourteen  characteristics  are 
considered  in  the  ULCAR  version  of  the  CHARS.FOR  program:  foF2, 
foFl,  Mfac,  MUF,  fmin,  foEs,  fmnF,  fmnE,  foE,  fxl,  hpF,  hpF2,  hpE, 
hpEs.  Definitions  of  these  characteristics  are  given  in  Tables  2  and  4. 

5.  Summary 

Programs  have  been  written  to  produce  the  monthly 
characteristic  database  files  for  uneven  time  sampling  as 
recommended  by  the  IIWG.  Supplemental  programs  have  been 
created  to  convert  the  original  ADEP  files  (type  1,  1988)  to  the 
standard  ADEP  output  files  (type  2,  1989)  and  from  these  to  the 
IIWG  ADEP  format  (type  3,  1990).  The  ULCAR  ADEP  program  now 
generates  files  in  the  IIWG  ADEP  format.  Any  monthly 
characteristics  file  in  this  format  can  be  input  to  the  characteristics 
program.  The  characteristics  program,  CHARS.FOR  is  given  in 
Appendix  A.  A  program  to  read  and  encode  the  characteristics 
database  file  is  given  in  Appendix  B.  This  program  can  be  easily 


modified  to  read  the  characteristics  data  for  use  in  studies  and  other 
applications. 

From  these  programs,  ionosonde  data  collected  from  the 
worldwide  Digisonde  network  [Reinisch  1986,  Tang  et.  al.  1990]  and 
from  other  networks  or  sounders  using  the  IIWG  ADEP  format  can 
quickly  and  confidently  be  used  to  generate  ionospheric 
characteristics  files  for  archiving.  The  later  use  of  these  files  is  easily 
accomplished  (see  Appendix  B). 
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APPENDIX  A 
PROGRAM  CHARS.FOR 


c.. .compile  with  fl/Gt3000 
$Storage;  2 
c 

PROGRAM  Chars 
C 

c... written  by  Bob  Gamache 

C  Read  ADEP  data  for  a  month  and  create  characteristics  files 
C  in  II WG  Workshop  format. 
c...July  10,  1989 

c...last  modified  December  13,  1990  by  Bob  Gamache 
c 

c...file  10  =  Standard  ADEP  output  for  one  month  of  data  (1989 
Format) 

c...file  3  =  IIWG  recommended  format  monthly  characteristics  data 

c...  (9000  byte  blocks) 

c...file  4  =  informative  data  for  the  month  -  check  feature 

c 

C 

parameter(N=31,  M=288) 
c 

Integer*2 

cfoF2,cfoFl,cM3000,cMUF,cfmin,cfoEs,cfminF,cfminE,cfoE, 

+  cfxI,chpF,chpF2,chpE,chpEs 
c 

character*! 

qfoF2,qMUF,qfoFl,qM3000,qfmin,qfoEs,qfminF,qfminE, 

+  qfoE,qfxI,qhpF,qhpF2,qhpE,qhpEs 
c 

Integer*4  Itimes.Ihr 
c 

INTEGER*2  IDn(80),IPREF(60),IAF(20),IOTF(400),IOThF(400), 

+ 

10  AF(400),IODF(400),IOTF1  ( 1 50),IOThFl  ( 1 50),lOAFl  ( 1 50),IODF1  (150 

) 

+ 

,IOTE(150),IOThE(150),IOAE(150),10DE(150),IXTF(400),lXAF(400), 

+ 

1XDF(400),IXTF1(150),IXAF1(150),IXDF1(150),1XTE(150),IXAE(150), 

+  IXDE(  1 50),MEDF(20),MEDE(20),MEDES(20),IedF(  1 20) 

REAL 

SC  ALED(45),GCONST(5),DTT(  1 6),FOF(400),FOF1  ( 1 50),FOE(  150), 

+  FXF(400),FXF1(150),FXE(150),THF2(20),THF1(20),THE(20), 

+  THMON(20),THVAL(20) 


Character*  120  Sysdes 
C 

CHARACTER *40  IFILEJout 
CHARACTER*30  Station,Systeni 
CHARACTER*  10  Scale(2),Edit(3) 

CHARACTER *5  cStat 
CHARACTER*  1 0  char(80),units(80),UT 
character*2  code(80) 
character*!  ql,aQL 
LOGICAL  EOF 

c... characteristics  data  arrays 

Common/Chard  1  /cfoF2(N,M),cMUF(N,M),cfoFl  (N,M) 
Common/chard2/cM3000(N,M),cfmin(N,M),cfoEs(N,M) 
Common/chard3/cfminF(N,M),cfminE(N,M),cfoE(N,M) 
Common/chard4/cfxI(N,M),chpF(N,M),chpF2(N,M) 
Common/chard5/chpE(N,M),chpEs(N,M) 
c 

c... characteristics  data  array  qualifiers 

Conimon/QLl/qfoF2(N,M),qMUF(N,M),qfoFl(N,M) 

Common/QL2/qM3000(N,M),qfmin(N,M),qfoEs(N,M) 

Common/QL3/qfmiinF(N,M),qfininE(N,M),qfoE(N,M) 

Common/QL4/qfxI(N,M),qhpF(N,M),qhpF2(N,M) 

Common/QL5/qhpE(N,M),qhpEs(N,M) 

c 

c... times  and  measurements 

Common/times/md(N),Ihv(N,24),NN,nb,Ndb(  1 5) 
Common/time4/Itimes(N,M) 
c 

dimension  aQL(0:l) 
c 

c...list  of  characteristics  followed  by  Units  and  URSI  codes 

data  char/*  foF2’,'  foFl','  M3000F2’,'  MUF3000F2’, 
+  '  fmin','  foEs','  fmnF’,'  fmnE’, 

+  '  foE','  fxr;  h"F,'  h’’F2', 

+  ’  h"E';  h"Es',  66*’  '/ 


c 


data  units/’  0.1  MHz’,’  0.01  MHz’,’  0.01’,’  0.1  MHz’, 

+  '0.1  MHz','  0.1  MHz','  0.01  MHz’,’  0.01  MHz’ 

+  ,'  0.01  MHz','  0.1  MHz’,’  1.0  km’,’  1.0  km’, 

+  ■  1.0  km',’  1.0  km’,  66*’  ’/ 


data  code/’OO’,’  1 0’,'03’,’07’,’42’,’30’,’80’,'8 1  ’,'20’,’5 1  ’,’  1 6’, 
’04',’24’,’34’,66*’  ’/ 
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c 

Scale(l)  =  '  Manual' 

Scale(2)  =  '  Automatic' 
c 

Edit(l)  =  '  Edited' 

Edit(2)  =  'Non-Edited' 

Edit(3)  =  '  Mixed' 
c 

c... qualifying  leters  from  ADEP  edit  flage 
aQL(l)  =  '  ' 
aQL(O)  =  'r 
lUNIT  =  10 
nl9  =  19 

c...#  characteristics 
rear  =  14 
c 

K  =  0 

C...DIGISONDE  operates  in  UT 
UT  =  '0000' 
c 

c.. .preset  arrays 
call  preset 
c 

WRITE(*,*)  'Enter  name  of  input  file.' 

READ(*,'(A)')  IFILE 

0PEN(UNIT=IUNIT,F1LE=IFILE,F0RM='F0RMATTED', 
+  ACCESS='SEQUENTIAL') 
c 

WRITE(*,*)  'Enter  name  of  output  file.' 

READ(*,'(A)')  lout 

OPEN(3,FILE=Iout,FORM='FORMATTED', 

+  ACCESS='SEQUENTIAL',  RECL=9000) 
c 

WRITE(*,*)  'Enter  ionosonde  station  name,  A30' 
read(*,'(a30)')  station 
c 

3  WRITE(*,*)  'Scaling  Type' 

Write(*,*)  '  -  Enter  1  for  Manual  scaling  ' 

Write(*,*)  '  -  Enter  2  for  Automatic  scaling  ' 

read(*,*)  Iscal 

if(Iscai  .ne.  1  .and.  Iscal  .ne.  2  )  then 
Write(*,*)  'incorrect  entry,  try  again' 
go  to  3 
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endif 


5  WRITE(*,*)  Data  Editing  ' 

Write(*,*)  '  -  Enter  1  for  Edited  ' 

Write(*,*)  '  -  Enter  2  for  Not-edited’ 

Write(*,*)  '  -  Enter  3  for  Mixed' 

read(*,*)  ledit 

if(Iedit  .ne.  1  .and.  ledit  .ne.  2  .and.  ledit  .ne.  3)  then 
Write(*,*)  'incorrect  entry,  try  again' 
go  to  5 
endif 
c 

WRITE(*,*)  'Enter  ionosonde  system  name,  A30' 
read(*,'(a30)')  System 
c 

ql  =  " 

if(Iedit  .eq.  2)  ql  =  7 
C 

15  CONTINUE 
K  =  K+  1 
CALL 

RDADEP(EOF,IUNIT,IDFI,IPREF,SCALED,GCONST,IAF,DTT,IOTF, 

+ 

I0ThF,I0AF,I0DF,F0F,I0TFl,I0ThFl,I0AFl,I0DFlJ70Fl,I0TE,I0ThE, 

+ 

10  AE,IODE,FOE,IXTF,IX  AF,IXDF,FXF,IXTF1  ,IXAF1  ,IXDF1  ,IXTE, 

+ 

IXAE,IXDE,FXE,MEDF,MEDE,MEDES,THF2,THFl,THE,THVAL,IedF,Sysdes) 

C 

if(EOF)  go  to  25 
c 

c...find  day  and  time 

lyr  =  IPREF(1)*10  +  IPREF(2) 
c...add  century  to  year 
lyr  =  n  19*  100  +  lyr 

Iday  =  IPREF(3)*100  +  IPREF(4)*10  +  IPREF(5) 

Ihr  =  IPREF(6)*10  +  IPREF(7) 

Min  =  IPREF(8)*10  +  IPREF(9) 

Isec  =  IPREF(10)*10  +  IPREF(ll) 
c...for  debuging 

write(l  l,'(5x,a,4i5)')  'ionogram  #,  yr,day,hour;',K,Iyr,Iday,Ihr 

c 

if(Iday.gt.366  .or.  Ihr.gt.23)  go  to  15 
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c 

c... month,  day  of  the  month,  and  number  of  days  in  the  month 
call  dmonth(lyr,  Iday,  Month,  Mday) 
if(Mday.gt.31)  go  to  15 
c 

c...for  hourly  values  the  day  starts  at  23:54 
if(Ihr  .eq.  23  .and.  Min  .gt.  54)  then 
Mday  =  Mday  +1 

c...do  not  exceed  days  in  the  month 
if(Mday  .gt.  NN)  go  to  15 
endif 
c 

c... increment  counter  for  number  of  measurements  per  day 
md(Mday)  =  md(Mday)  +  1 
if(md(Mday)  .gt.  288) 

+  write(*,*)  ’  Warning  md(',Mday,’)=)',md(Mday) 
c...time  of  measurement 

Itimes(Mday,md(mday))  =  Ihr*  10000  +  Min*  100  +  Isec 
c 

c... store  characteristic  data  in  arrays 

cfoF2(Mday,md(Mday))  =  ((scaled(l)+.001)*10) 
qfoF2(Mday,md(Mday))  =  aQL(IedF(l)) 
cfoFl(Mday,md(Mday))  =  ((scaled(2)+.001)*10) 
qfoFl(Mday,md(Mday))  =  aQL(IedF(2)) 
if(scaled(24).eq.3000.)  then 
if(scaled(3)  .It.  50.)  then 

cM3000(Mday,md(Mday))  =  ((scaled(3)+.0001)*100) 

else 

cM3000(Mday,md(Mday))  =  9999 
endif 

cMUF(Mday,md(Mday))  =  ((scaled(4)+.001)*10) 

else 

cM3000(Mday,md(Mday))  =  9999 
cMUF(Mday,md(Mday))  =  9999 
endif 

qM3000(Mday,md(Mday))  =  aQL(IedF(3)) 
qMUF(Mday,md(Mday))  =  aQL(IedF(4)) 
cfmin(Mday,md(Mday))  =  ((scaled(5)+.001)*10) 
qfmin(Mday,md(Mday))  =  aQL(IedF(5)) 
cfoEs(Mday,md(Mday))  =  ((scaled(6)+.001)*10) 
qfoEs(Mday,md(Mday))  =  aQL(IedF(6)) 
cfminF(Mday,md(Mday))  =  ((scaled(7)+.001)*10) 
qfminF(Mday,md(Mday))  =  aQL(IedF(7)) 
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cfminE(Mday,md(Mday))  =  ((scaled(8)+.001)*10) 
qfminE(Mday,md(Mday))  =  aQL(IedF(8)) 
cfoE(Mday,md(Mday))  =  ((scaled(9)+.001)*10) 
qfoE(Mday,md(Mday))  =  aQL(IedF(9)) 
cfxI(Mday,md(Mday))  =  ((scaled(10)+.001)*10) 
qfxI(Mday,md(Mday))  =  aQL(IedF(10)) 
chpF(Mday,md(Mday))  =  ((scaled(l  1)+.01)*1) 
qhpF(Mday,md(Mday))  =  aQL(IedF(ll)) 
chpF2(Mday,md(Mday))  =  ((scaled(12)+.01)*l) 
qhpF2(Mday,md(Mday))  =  aQL(IedF(12)) 
chpE(Mday,md(Mday))  =  ((scaled(13)+.01)*l) 
qhpE(Mday,md(Mday))  =  aQL(IedF(13)) 
chpEs(Mday,md(Mday))  =  ((scaled(14)+.01)*l) 
qhpEs(Mday,md(Mday))  =  aQL(IedF(14)) 

go  to  15 

25  continue 

..station  number,  latitude  N,  longitude  E 
iScode  =  int(GCONST(5)) 

Rlat  =  GCONST(3) 

Rlong  =  GCONST(4) 

Write(cStat,'(2x,i3) )  iScode 

..calculate  hourly  values  for  medians 
call  hourly 

write(4,*)  'measurements/day' 
write(4,'(31I4,//)')  md 

..find  total  number  of  measurements 
nmt  =  0 
do  30  L=1,NN 
30  nmt  =  nmt  +  md(L) 

..calculate  the  number  of  day.s  to  write  into  a  block 
call  blocker 

..output  first  block 

write(3,’(A30,A5,A4,2F5.1,2A10,A30)’)  station,cStat,UT,Rlat, 
+  Rlong, Scale(Iscal),Edit(Iedit),System 
write(3,'(30I4)')  Iyr,Month,NN,ncar,nmt,md 
write(3,'(12A10y)  (Char(L),L=l,ncar) 


write(3,'(12A10)’)  (units(L),L=l,ncar) 
write(3,’(60A2)’)  (code(L),L=l,ncar) 

.output  measurement  times. .depends  on  the  number  of  blocks 
call  outime 

.compute  mendian  points 

call  findMed(cfoF2,  qfoF2,  Medql,  900) 
call  output(cfoF2,  qfoF2,  Medql,  900) 

call  findMed(cfoFl,  qfoFl,  Medql,  900) 
call  output(cfoFl ,  qfoFl,  Medql,  900) 

call  findMed(cM3000,  qM3000,  Medql,  900) 
call  output(cM3000,  qM3000,  Medql,  900) 

call  findMed(cMUF,  qMUF,  Medql,  900) 
call  output(cMUF,  qMUF,  Medql,  900) 

call  findMed(cfmin,  qfmin,  Medql,  900) 
call  output(cfmin,  qfmin,  Medql,  900) 

call  findMed(cfoEs,  qfoEs,  Medql,  900) 
call  output(cfoEs,  qfoEs,  Medql,  900) 

call  findMed(cfminf,  qfminf,  Medql,  900) 
call  output(cfminf,  qfminf,  Medql,  900) 

call  findMed(cfminE,  qfminE,  Medql,  900) 
call  output(cfminE,  qfminE,  Medql,  900) 

call  findMed(cfoE,  qfoE,  Medql,  900) 
call  output(cfoE,  qfoE,  Medql,  900) 

call  findMed(cfxI,  qfxl,  Medql,  900) 
call  output(cfxI,  qfxl,  Medql,  900) 

call  findMed(chpF,  qhpF,  Medql,  900) 
call  output(chpF,  qhpF,  Medql,  900) 

call  findMed(chpF2,  qhpF2,  Medql,  900) 
call  output(chpF2,  qhpF2,  Medql,  900) 


call  findMed(chpE,  qhpE,  Medql,  900) 
call  output(chpE,  qhpE,  Medql,  900) 
c 

call  findMed(chpEs,  qhpEs,  Medql,  900) 
call  output(chpEs,  qhpEs,  Medql,  900) 
c 

END 

c 

subroutine  output(Icharx,  qldl,  Medql,  Ilimit) 
c 

parameter(N=3 1 ,  M=288) 
c 

character*5  A5out(24) 
character*70  A70 
character*!  qldl,qq,qa 
character*!  ql,  dl 
c 

c.. .times  and  measurements 

Common/times/md(N),Ihv(N,24),NN,nb,Ndb(  1 5) 
c 

common/median/Med(24),  Icount(24),  nUQ(24),  LQ(24),  nUD(24), 
+  LD(24),  IRange(24) 
c 

dimension  Icharx(N,M),  qldl(N,M) 
c 

qa  =  7 

c... output  characteristic,  medians,  counts,  upper/lower  quartiles, 
c.. .deciles  and  range 
c 

c...set  qualifying  and  descriptive  letters  for  medians 
if(Medql.eq.l)  then 
ql  =  7 
dl  =  " 
else 
ql  =  7 
dl  =  7 
endif 

c... output  characteristics 
c 

call  blank(A5out) 
ns  =  0 

c...loop  number  of  days 
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do  10  nd  =  1,  NN 
c...loop  measuremnets  per  day 
do  9  I  =  1,  md(nd) 
ns  =  ns  +  1 
c 

if(Icharx(nd,I)  .gt.  Ilimit  .or. 

+  Icharx(nd,I)  .eq.  0)  then 

A5out(ns)  =  ' 
else 

qq  =  qldl(nd,I) 

write(A5out(ns),'(I3,2Al)’)  Icharx(nd,l),qq,qq 

endif 

c...when  full  output  the  start  again 
if(ns  .eq.  24)  then 

WTite(3,'(24A5)')  ( A5out(kk),kk=l  ,ns) 
ns  =  0 
endif 
c 

9  continue 

10  continue 
c 

c... write  out  if  there  is  any  remainding  data 

if(ns  .gt.  0)  write(3,’(24A5)’)  (A5out(kk),kk=l,ns) 
c 

c... medians 

do  15  L  =  1,  24 

write(A5out(L),’(I3,2Al )')  Med(L),ql,dl 
15  continue 

write(3,'(24A5)’)  (A5out(L),L=l  ,24) 
c 

c.. .median  counts  and  Range 
do  20  L  =  1,  24 

write(A5out(L),'(I2,I3)')  Icount(L),IRange(L) 

20  continue 

write(3,'(24A5)')  (A5out(L),L=l  .24) 
c 

c... Upper  Quartiles 
do  25  L  =  1,  24 

write(A5out(L),’(I3,2AlV)  nUQ(L),ql,dl 
25  continue 

write(3,’(24A5)')  ( A5out(L),L=l ,24) 
c 

c... Lower  Quartiles 
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do  30  L  =  1,  24 

write(A5out(L),’(13,2Al)*)  LQ(L),ql,dl 
30  continue 

write(3,’(24A5)’)  ( A5out(L).L=l ,24) 
c 

c.. .Upper  Deciles 

do  35  L  =  1,  24 

writeC A5out(L),’(I3 .2 A 1  )*)  nUD(L),ql,dl 
35  continue 

write(3 ,’(24 A5 )')  ( A5out(L),L=l  ,24) 
c 

c... Lower  Deciles 

do  40  L  =  1,  24 

write(A5out(L),’(I3,2Al)*)  LD(L),ql,dl 
40  continue 

write(3,’(24A5)’)  (A5out(L),L=l  ,24) 
c 

50  continue 
return 
end 
c 

subroutine  outime 
c 

parameter(N=3 1 ,  M=288) 
c 

character*6  A6out(20) 
c 

c... times  and  measurements 
Integer*4  Itimes 
Common/time4/Itimes(N,M) 
Common/times/md(N),Ihv(N,24),NN,nb,Ndb(  15) 
c 

c.. .output  measurement  times 
c 

call  blank(A6out) 
ns  =  0 

c...loop  number  of  days 
do  10  nd  =  1,  NN 
c...loop  measuremnets  per  day 
do  10  I  =  1,  md(nd) 
ns  =  ns  +  1 

write(A6out(ns),'(I6)')  Itimes(nd,I) 
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c 


if(ns  .eq.  20)then 

write(3,'(20A6)')  ( A6out(kk),kk=  1  ,ns) 
ns  =  0 
endif 

10  continue 

c 

c... write  out  if  there  is  any  remainding  data 

if(ns  .gt.  0)  write(3,'(20A6)’)  (A6out(kk),kk=l,ns) 
c 

return 

end 

c 

subroutine  hourly 

c 

parameter(N=3 1 ,  M=288) 
c 

Integer*4  Itimes,Ihr,min 
c 

c... times  and  measurements 

Common/times/md(N),Ihv(N,24),NN,nb,Ndb(  1 5) 
Common/time4/Itimes(N,M) 
dimension  Jdel(N,M) 
c 

do  5  1=1,31 
do  5  J=l,24 
5  JdeKU)  =  10 
c 

c...set  hourly  times 
c...loop  days  -  I 
do  50  1=1,  NN 

c...loop  measurements/day  --  J 
do  49  J=l,  md(I) 

Jmk  =  J 
c 

Ihr  =  Itimes(l,J)/10000 
Min  =  (Itimes(I,Jy  -  Ihr*  10000)/! 00 
c 

if(Min  .gt.  54)  then 
Mk  =  Ihr  +  1 
Id  =  I 
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Idel  =  60  -  Min 
go  to  20 
endif 
c 

if(Min  .It.  6)  then 
Mk  =  Ihr 

if(Mk.eq.O)  Mk  =  24 
Id  =  I 
Idel  =  Min 
go  to  20 
endif 
c 

c...not  in  proper  time  zone 
go  to  49 
c 

c... hourly  value  check  for  closest  value 
20  continue 

if(Ihv(Id,Mk)  .eq.  0)  then 
c...set  initial  value 

Ihv(Id,Mk)  =  Jmk 
Jdel(Id.Mk)  =  Idel 
else 

c. ..compare  to  previous  value 

if(Idel  .le.  Jdel(Id,Mk))  then 
Ihv(Id,Mk)  =  Jmk 
Jdel(Id,Mk)  =  Idel 
endif 
endif 
c 

49  continue 

50  continue 
c 

write(4,'(/,2x.24I3)')  ((Ihv(I J).J=1 ,24)J=1 ,3 1 ) 
return 
end 
c 

^  Kc  Hr  Dr  4ii|i  #  4i  4i  4i  4i «  Hi  Hi  *  4i  4c  4i  4i 

subroutine  Blank(AA) 
c 

character*  120  A  A 
ns  =  1 

do  10  I  =  1,  12 
AA(ns:ns+9)=’ 
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10  ns  =  ns  +  10 
return 
end 


c 

subroutine  findMed(Icharx,  qcharx,  Medql,  Ilimit) 
c 

parameter(N=3 1 ,  M=288) 
c 

c... times  and  measurements 

Common/times/md(N),Ihv(N,24),NN,nb,Ndb(  1 5) 


c...find  medians,  counts,  upper/lower  quartiles-deciles  and  range 
character*!  qcharx 
c 

common/median/Med(24),  Icount(24),  nUQ(24),  LQ(24),  nUD(24), 
+  LD(24),  IRange(24) 

dimension  Icharx(N,M),  qcharx(N,M),  Ival(N), 

+  Ic  1  (3 1  ),lc2(3 1  ),mc  1(31  ),mc2(3 1 ) 


c...lcl,  lc2,  mcl,  and  mc2  are  the  URSI  location  vectors  for  the  upper 
c...and  lower  quartiles  (see  UAG-23  pl85) 

data  lcl/1 , 1 , 1 , 1 , 1 ,2,2,2,2,3,3,3,3,4,4,4,4,5,5,5,5,6,6,6,6,7,7,7, 

+  7,8,8/ 

data  lc2/l,l,2,2,2,2,2,3,3,3,3,4,4,4,4,5,5,5,5,6,6,6,6,7,7,7,7,8, 

+  8,8,8/ 

data  mcl/1, 1,2,3,4,5,6,6,7,8,9,9,10,11, 12,12,13, 14,15,15, 16,17, 

+  18,18,19,20,21,21,22,23,24/ 

data  mc2/l,l,3,4,5,5,6,7,8,9,9,10,l  1,1 1,12,13,14,14,15,16,17,17, 
+  18,19,20,20,21,22,23,23,24/ 


c 

c... preset 

do  2  1=1,24 
Med(I)  =  0 
Icount(I)  =  0 
Irange(I)  =  0 
nUQ(I)  =  0 
LQ(I)  =  0 
nUd(I)  =  0 
LD(1)  =  0 
2  continue 
Medql  =  1 
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c 


c...loop  over  hours 
do  50  I  =  1,  24 
c 

do  3  J=l,31 
Ival(J)  =  0 
3  continue 
c 

c...get  hourly  value  for  each  day 
nc  =  0 
c 

do  45  J  =  1,  NN 
Hoc  =  Ihv(J,I) 
if(Iloc  .ne.  0)  then 

if(Icharx(J,Iloc)  .It.  Ilimit  .and. 

+  Icharx(J,Iloc)  .ne.  0)  then 

nc  =  nc  +  1 

Ival(nc)  =  Icharx(J,Iloc) 
if(qcharx(J,Iloc)  .eq.  7)  Medql  =  0 

endif 

endif 

45  continue 
c 
c 

c...set  counter 

Icount(I)  =  nc 
c 

if(nc  .eq.  0)  go  to  50 
c.. .order  the  values 

call  Order(Ival,  nc) 
c 

c...nc=l 

if(nc  .eq.  1)  then 
Med(I)  =  Ival(l) 
go  to  50 
endif 
c 

c...nc>2 

if(nc  .ge.  2)  then 
c 

c... median  is  at  nc/2 
nc2  =  nc/2 

if(Mod(nc,2)  .eq.  0)  then 
Med(i)  =  (Ival(nc2)  +  Ival(nc2+1))  /  2 
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else 

Med(I)  =  Ival(nc2+1) 
endif 
c 

c...quartiles 

if(nc  .ge.  3)  then 

nUQ(I)  =  (Ival(mcl(nc))  +  Ival(mc2(nc))  )/2 
LQ(I)  =  (Ival(lcl(nc))  +  Ival(lc2(nc))  )/2 
IRange(I)  =  nUQ(I)  -  LQ(I) 
endif 
c 

c.. .deciles 

if(nc  .ge.  10)  then 
nclO  =  nc/10 
nUD(I)  =  IvaI(nc-nclO+l) 

LD(I)  =  Ival(nclO) 
endif 
c 

c...nc>2 

endif 

c 

50  continue 
return 
end 
c 

subroutine  Order(Kval,  nc) 
c 

c...sort  integer  array  Kval  in  assending  order 
dimension  Kval(31),  Idum(31) 
c 

c... there  are  nc  values 
c 

Kount  =  nc 
do  100  I  =  1,  nc 
max  =  0 
ip  =  0 

c...find  max  value 
do  50  J  =  1 ,  nc 
if(Kval(J)  .ge.  Max)  then 
ip  =  J 

Max  =  Kval(J) 
endif 
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50  continue 

Idum(Kount)  =  Max 
Kval(ip)  =  0 
Kount  =  Kount  -  1 
100  continue 
c 

c...put  ordered  array  back  into  Kval 
do  150  K  =  I,  nc 
150  Kval(K)  =  Idum(K) 
c 

return 

end 

c 

subroutine  dmonth(Iyr,  Iday,  Month,  Mday) 
c 

parameter(N=3 1 ,  M=288) 
c 

c... times  and  measurements 

Common/times/md(N),Ihv(N,24),NN,nb,Ndb(15) 

c 

c...from  Iday  calculate  month,  day  of  the  month,  and  number 
c...of  days  in  the  month, 
c 

c...Mth  first  day  of  each  month  (Feb  start) 

c...Leap  leap  year  correction  factor 

c...Iyr  year  (12)  e.g.  83  corresponds  to  1983 

c.. .Month  local  month 

c...Iday  local  day  (1-365) 

c...Mday  local  day  (1-31) 

c...Nmonth  #  of  days  for  each  month 

C...NN  #  of  days  in  the  current  month 

c 

dimension  Mth(12),Nmonth(12) 
c 

data  Mth/32,60,91, 121, 152,182,213,244,274,305,335,366/ 
data  Nmonth/3 1,28, 3 1,30,3 1,30,3 1,3 1,30,  31,  30,  31/ 
c 

c...leap  year  correction 
Leap  =  0 

if(MOD(Iyr,4)  .eq.  0)  Leap  =  I 
c 

if(Iday  .It.  32  )  then 
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I  =  1 

Month  =  1 
Mday  =  Iday 
NN  =  Nmonth(Month) 
go  to  999 
endif 
c 

do  10  I  =  2,  12 

if(Iday  .It.  Mth(I)+Leap)  go  to  11 

10  continue 
c 

11  Month  =  I 
if(Month  .gt.  2)  then 

Mday  =  Iday  -  Mth(Month-l)  -  leap  +  1 
NN  =  Nmonth(Month) 
else 

Mday  =  Iday  -  Mth(Month-l)  +  1 
NN  -  Nnionth(Month)  +  Leap 
endif 
c 
c 

999  return 
end 
c 

subroutine  blocker 
c 

parameter(N=31,  M=288) 


c 

c...md(i) 

number  of  measurements  on  day  i 

C...NN 

number  of  days  in  the  currenbt  month 

c...nb 

number  of  blocks  per  characteristic 

c...ibc 

max  number  of  characters  per  block 

c...nc 

number  of  characters  per  block 

c...Idb 

counter  for  number  of  days/block 

c...Ndb(15)  number  of  days  in  each  block  to  be  written,  15  would 
c...  support  a  2  minute  measurement  schedule 

c 

c... times  and  measurements 

Common/times/md(N),Ihv(N,24),NN,nb,Ndb(15) 

c 

ibc  =  9000 
nc  =  0 
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c.. .preset 
nb  =  1 
Idb  =  0 
do  5  i=l,15 
5  Ndb(i)  =  0 
Ndb(l)  =  31 
c 

do  20  i=l,NN 

c...six  characters  per  measurement 
Idb  =  Idb  +  1 
nc  =  nc  +  md(i)*6 
if(nc  .gt.  ibc)  then 
Ndb(nb)  =  Idb  -  1 
Idb  =  1 
nb=nb  +  1 
nc  =  md(i)*6 
endif 

20  continue 

c...when  more  than  one  block  remainder  goes  in  last  block 
if(nb  .gt.  1)  Ndb(nb)  =  Idb 
c 

write(4,*)  'nb,  Ndb' 
write(4,'(I4,5x,15I4,//)’)  nb,  Ndb 
return 
end 
c 

^  ,tc  ^  :|c  :<(  ik  Ik  *  i)»l<  3f<  *  i)<  *  *  Ik 

subroutine  preset 
c 

parameter(N=3 1 ,  M=288) 
c 

Integer*2 

cfoF2,cfoFl,cM3000,cMUF,cfmin,cfoEs,cfminF,cfminE,cfoE, 

+  cfxI,chpF,chpF2,chpE,chpEs 
c 

Common/Chard  1  /cfoF2(N,M),cMUF(N,M),cfoFl  (N,M) 
Common/chard2/cM3000(N,M),cfmin(N,M),cfoEs(N,M) 
Common/chard3/c^minF(N,M),cfminE(N,M),cfoE(N,M) 
Common/chard4/cfxI(N,M),chpF(N,M),chpF2(N,M) 
Common/chard5/chpE(N,M),chpEs(N,M) 
c 

c... times  and  measurements 

Common/times  1  /md(N),Ihv(N,24),NN,nb,Ndb(  1 5) 
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Integer*4  Itimes 
Coiiimon/time4Atimes(N,M) 
c 

do  3  k=l,24 
do3  i=l,N 
Ihv(i,k)  =  0 
3  md(i)  =  0 
c 

do  10  i=l,N 
do  10  j=l,M 
Itimes(i,j)  =  0 
cfoF2(i  j)  =  9999 
cfoFl(io)  =  9999 
cM3000(i,j)  =  9999 
cMUF(io)  =  9999 
cfmin(i,j)  =  9999 
cfoEs(iJ)  =  9999 
cfminF(i,j)  =  9999 
cfniinE(i,j)  =  9999 
cfoE(i,j)  =  9999 
cfxl(i,j)  =  9999 
chpF(i,j)  =  9999 
chpF2(i,j)  =  9999 
chpE(i,j)  =  9999 
10  chpEs(i,j)  =  9999 
c 

return 

end 

C 

c...read  1990  ADEP  IIWG  suggested  format 
Subroutine 

RDADEP(EOF,IU,IDn,IPREF,SCALED,GCONST,IAF,DTT,IOTF, 

+ 

I0ThF,10AF,I0DF,F0F,I0TFl  JOThFl  ,IOAFl  ,IODFl  ,FOFl  ,10TE,10ThE, 

+ 

lOAEJODE  J?OE,IXTF,IXAF,IXDF,FXFJXTFl  ,1X  AFl  JXDFl  ,FXF1  ,IXTE, 

+ 

IXAE,IXDE,FXE,MEDF,MEDE,MEDES,THF2,THFl,THE,THVAL,IedF,Sysdes) 

C 

C  This  subroutine  reads  scaled  Ionospheric  data  from  a  textfile 

C  in  the  format  specified  by  the  University  Of  Lowell  Center  for 
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C  Atmospheric  Research  for  ionospheric  data  produced  by 
Digisonde 

C  256  digital  ionospheric  sounders  and  processed  by  the  ADEP 
C  ARTIST  Data  Editing  Program.  The  format  is  that  suggested  by 
the 

C  IIWG  Workshop  --  ADEP  90  - 

C  The  EOF  variable  is  set  if  an  End  Of  File  is  encountered  during 
C  the  file  read. 

C 

C  For  further  information,  see  the  ULowell  doccument: 

C  ARTIST  Data  Editing  Program  Output  Format  (02Apr88) 

C  NOTE:  Values  which  are  all  conatined  on  one  line  are  not 
C  read  in  implied  DO  loops  so  as  to  minimize  the  effect 

C  of  errors  on  the  subsiquent  records. 

C 

C  Reads  data  from  the  FORTRAN  unit  number  lU  into  the  data 
arrays. 

C 

LOGICAL  EOF 
INTEGER  I,IU 
INTEGER 

IDFI(80),IPREF(60),IAF(20),IOTF(400),IOThF(400),IOAF(400), 

+ 

IODF(400),IOTF1  ( 1 50),IOThFl  ( 1 50),IOAF1  ( 1 50).IODF1  ( 1 50),IOTE(1 50 

) 

+ 

,IOThE(  1 50),IO  AE(  1 50),IODE(  1 50),IXTF(400),IX  AF(400),IXDF(400), 

+  IXTF1(150),IXAFI(150),IXDF1(150),1XTE(150),IXAE(150), 

+  IXDE(150),MEDF(20),MEDE(20),MEDES(20),IedF(120) 

REAL 

SCALED(30),GCONST(5^,DTT(16),FOF(400),FOF1(150),FOE(150), 

+  FXF(400),FXF1  ( 1 50),FXE(  1 50),THF2(20),THF1  (20),THE(20), 

+  THMON(20),THVAL(20) 

Character*  120  Sysdes 
C 

c... formats 

101  FORMAT  (4013) 

102  FORMAT  (16F7.3) 

103  Format(A120) 

104  FORMAT  (120Z1) 

105  FORMAT  (15F8.3) 

106  FORMAT  (6012) 

107  FORMAT  (12011) 
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108  FORMAT  (20F6.3) 

109  FORMAT  (13E9.4E1) 
c...data  file  index 

C  The  data  file  index  integers  should  all  be  on  one  line. 
READ(IU,101,END=9)  IDH 

..Geophysical  constants  --  Code  1 
IF(IDFI(1).GT.0)  READ(IU,102,END=9)  GCONST 

..System  description  —  Code  2 

if(IDFI(2).gt.O)  READ(IU,103,END=9)  Sysdes 

...ionogram  sounding  settings  (preface)  —  Code  3 

if(IDn(3).gt.O)  READ(IU,104.END=9)  (IPREF(I),I=1,IDFI(3)) 

...scaled  ionogram  parameters 

IF(IDFI(4).GT.O)  READ(IU,105,END=9)  (SCALED(I),I=1,IDFI(4)) 

ARTIST  analysis  flags 
IF(IDn(5)  .GT.O)  READ(IU,106,END=9)  IAF 

...Doppler  translation  table 

IF(IDn(6).GT.O)  READ(IU,102,END=9)  DTT 

...0-trace  F2  points 
..virtual  height 

IF(IDFI(7).GT.O)  READ(IU,101,END=9)  (IOTF(I),I=l,IDFI(7)) 
..true  height 

IF(IDFI(8).GT.O)  READ(IU,101,END=9)  (IOThF(I),I=l,IDFI(8)) 
..Amplitudes 

IF(IDFI(9).GT.O)  READ(IU,106,END=9)  (IOAF(I),I=l,IDFI(9)) 
..Doppler  numbers 

IF(IDFI(10).GT.0)  READ(IU,107,END=9)  (IODF(I),I=1,IDFI(10)) 
..Frequency  table 

IF(IDFI(11).GT.0)  READ(IU,108,END=9)  (FOF(I),  I=l,IDn(ll)) 

...0-trace  FI  points 
c.. .virtual  height 

IF(IDFI(12).GT.O)  READ(IU,101,END=9)  (IOTFl(I),I=l,IDFI(12)) 
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C 

c...true  height 

IF(IDFI(13).GT.O)  READ(IU,101,END=9)  (I0ThFl(l),I=l,IDn(13)) 
c 

c... Amplitudes 

IF(IDFI(14).GT.O)  READ(IUJ06,END=9)  (I0AF1(I),I=1,IDFI(14)) 
c 

c... Doppler  number 

IF(IDFI(15).GT.O)  READ(IU,107,END=9)  (I0DF1(I),I=1,IDFI(15)) 
c 

c.. .Frequency  table 

IF(IDFI(16).GT.O)  READ(IU,108,END=9)  (FOFl(I),  I=l,IDn(16)) 

...0-trace  E  points 
..virtual  heights 

IF(IDn(17).GT.O)  READ(IU,101,END=9)  (I0TE(I),I=1,IDFI(17)) 
..true  height 

IF(IDFI(18).GT.O)  READ(IU,101,END=9)  (I0ThE(I),I=l,IDFI(18)) 
..Amplitudes 

IF(IDFI(19).GT.O)  READ(IU,106,END=9)  (I0AE(I)J=1,IDFI(19)) 
..Doppler  numbers 

IF(IDFI(20).GT.0)  READ(IU,107,END=9)  (IODE(I),I=1,IDFI(20)) 
..Frequency  table 

IF(iDFI(21).GT.O)  READ(IU,108,END=9)  (FOE(I),  I=1,IDFI(21)) 

...X -trace  F2  points 
..virtual  heights 

IF(IDFI(22).GT.O)  READ(IU,101,END=9)  (IXTF(I),I=1,IDFI(22)) 
..Amplitudes 

IF(rDFI(23).GT.O)  READ(IU,106,END=9)  (IXAF(I),I=1,IDFI(23)) 
..Doppler  numbers 

IF(IDn(24).GT.O)  READ(IU,107,END=9)  (IXDF(I),I=l,IDn(24)) 
..Frequency  table 

IF(IDFI(25).GT.O)  READ(IU,108,END=9)  (FXF(I),I=1,IDFI(25)) 
..X-trace  FI  points 
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c... virtual  heights 

IF(IDFI(26).GT.O)  READ(IU,101,END=9)  (IXTF1(I),I=1,IDFI(26)) 
c 

c... Amplitudes 

IF(IDFI(27).GT.O)  READ(IU,106,END=9)  (1XAF1(I),1=1,IDF1(27)) 
c 

c... Doppler  numbers 

IF(IDFI(28).GT.O)  READ(IU,107,END=9)  (IXDF1(I),I=1,IDFI(28)) 
c 

c.. .Frequency  table 

IF(IDFI(29).GT.O)  READ(IU,108,END=9)  (FXF1(I),1=1,IDFI(29)) 

...X-trace  E  points 
..virtual  heights 

IF(IDFI(30).GT.0)  READ(IU,10l,END=9)  (IXTE(I),I=1,IDFI(30)) 
..Amplitudes 

IF(IDFI(31).GT.O)  READ(IU,106,END=9)  (IXAE(I),I=1,IDF1(31)) 
..Doppler  numbers 

IF(IDFI(32).GT.O)  READ(IU,107,END=9)  (IXDE(I),I=1,IDFI(32)) 
..Frequency  table 

IF(IDFI(33).GT.O)  READ(IU,108,END=9)  (FXE(I),I=1,IDFI(33)) 
...Median  amplitude  of  F  echo 

IF(IDn(34).GT.O)  READ(IU,106,END=9)  (MEDF(I),I=1  JDFI(34)) 
...Median  amplitude  of  E  echo 

IF(IDFI(35).GT.O)  READ(IU,106.END=9)  (MEDE(I),I=1,IDFI(35)) 
...Median  amplitude  of  Es  echo 

IF(IDFI(36).GT.O)  READ(IU,106,END=9)  (MEDES(I),I=1,IDFI(36)) 

...F2  layer  true  height  parameters 

IF(IDFI(37).GT.O)  READ(IU,109,END=9)  (THF2(I),I=1,IDFI(37)) 

...FI  layer  true  height  parameters 

IF(IDFI(38).GT.O)  READ(IU,109,END=9)  (THF1(I),I=1,IDFI(38)) 
...E  layer  true  height  parameters 

IF(IDFI(39).GT.O)  READ(IU,109,END=9)  (THE(I),I=1,IDFI(39)) 

...Valley  parameters 

IF(IDFI(40).GT.0)  READ(IU,109,END=9)  (THVAL(I),I=1,IDF1(40)) 
C 


c...Edit  Flags 

IF(IDFI(41).GT.O)  READ(IU,107.END=9)  (IedF(I),I=l,IDFI(41)) 
c 

EOF  =  .FALSE. 

RETURN 

9  EOF  =  .TRUE. 

RETURN 

END 

c 
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APPENDIX  B 

PROGRAM  RD-CHARS.FOR 
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Storage;  2 
c 

PROGRAM  Rd_Chars 
C 

c... written  by  Bob  Gamache 
C  Reads  database  characteristics  files 
C  in  IIWG  Workshop  format. 

..last  modified  13  December,1990  by  Bob  Gamache 

parameter(N=3 1 ,  M=288) 
c 

CHARACTER*40  IFILE.Iout 
CHARACTER*30  Station.System 
CHARACTER*  10  Scale,Edit 
CHARACTER*5  cStat,Ch5(9000).C5(24) 

CHARACTER*6  Ch6(9000) 

CHARACTER*  10  char(80),units(80),UT 
character*2  code(80) 
dimension  Icount(24),IRange(24) 
c 

c... times  and  measurements 

Common/tm/md(N),NN,nmt 

c 

WRITE(*,*)  '  Enter  name  of  input  file.’ 

READ(*;(A)’)  IFILE 

OPEN(  1  ,FILE=IFILE,FORM='FORMATTED’, 

+  ACCESS='SEQUENTIAL',RECL=9000) 
c 

WRITE(*,*)  'Enter  name  of  output  file.' 

READ(*,’(A)')  lout 

OPEN(2,FILE=Iout,FORM=’FORMATTED’, 

+  ACCESS='SEQUENTIAL',  RECL=9000) 
c 

c... read/write  first  record 

read(  1 ,'( A30, A5,A4,2F5. 1 ,2A  1 0,A30)',end=500)  station,cStat,UT, 
+  Rlat,Rlong,Scale,Edit,System 

read(  1  ,'(3014)’)  Iyr,Month,NN,ncar,nmt,md 
read(  1 ,'( 1 2A 1 0)')  (Char(L),L=  1  .near) 
read(  1 ,’( 1 2A 10)')  (units(L),L=l  ,ncar) 
read(  1  ,'(60A2)')  (code(L),L=l  ,ncar) 
c 
c 

write(2,'(A30,5x,A5,a,f5.1,a,f5.1,5x,a,A4)')  Station, cStat, 
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+’  latitude:', Rlat,'  longitude:',Rlong,'  UT:',UT 

write(2,'(lx,al0,5x,al0,2x,a,A30,a,I2,lx,i4)’)Scale,Edit, 

+'  System:  ’.System,'  Month:  ',Month,Iyr 
c 

Write(2,’(/,I5,a,I5,a)')  NN,'  Days  ',nmt,’  #  measurements  total' 
Write(2,'(/,a)')  '  #  measurements  per  day' 

Write(2,’(30I4)’)  md 
c 
c 

write(2,'(/,a,/)’)'  measurement  times:  hours:minutes:seconds' 
c... output  measurement  times 
c.. .number  of  measurement  lines 
nml  =  nmt/20 

if(  (nmt-nml*20)  .gt.  0)  nml  =  nml  +  1 
c... number  in  last  line 
nil  =  nmt  -  nml*20 
c...read  20  words  at  one  time 
ns  =  1 

ne  =  ns  +  19 
c 

do  20  L=l,  nml 

read(  1  ,'(20A6)’)  (Ch6(lck),kk=ns,ne) 
ns  =  ns  +  20 
ne  =  ns  +  19 
c...last  line 

if(l  ,eq.  nml)  ne  =  ns  +  nil  -  1 
20  continue 

write(2,'(15(2x,A6))')  (CH6(j),  j=l,  nmt) 
c 

c... output  characteristics 
c 

do  40  Nc=l,  near 

write(2,'(/,a,al0,lx,al0,2x,al  l,a2/)')’  characteristic:  ', 

+  char(Nc),units(Nc),'  URSl  code  ',code(Nc) 
c 

c.. .number  of  measurement  lines  for  characteristics 
nml  =  nmt/24 

if(  (nmt-nml*24)  .gt.  0)  nml  =  nml  +  1 
c... number  in  last  line 
nil  =  nmt  -  nml*24 
c...read  24  words  at  one  time 
ns  =  1 

ne  =  ns  +  23 
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c 


do  30  L=l,  nml 

read(  I  ,'(24A5)’)  (Ch5(kk),kk=ns,ne) 
ns  =  ns  +  24 
ne  =  ns  +  23 
c...last  line 

if(l  .eq.  nml)  ne  =  ns  +  nil  -  1 
30  continue 

write(2,'(15(3x,A5))’)  (ChSO),  j=l.  nmt) 
c 

c.. .medians 

read(l,'(24A5)')  (C5(L),L=1,24) 
write(2,'(/,a,/)')'  medians' 
write(2,'(12(3x,A5))')  (C5(J),  J=l,  24) 
c 

c... median  counts  and  Range 

read(l, ’(24(12, 13))')  (Icount(L),IRange(L),L=l,24) 
write(2,’(/,a,/)')'  counts’ 

write(2,’(12(3x,I2))’)  (Icount(J),  J=l,  24) 
c 

c... Range 

write(2,’(/,a,/)')'  range’ 

write(2,’(12(3x,I3))’)  (Irange(J),  J=l,  24) 
c 

c... Upper  Quartiles 

read(l,’(24A5)')  (C5(L),L=1,24) 
write(2,’(/,a,/)')'  upper  quartiles’ 

write(2,’(12(3x,A5))’)  (C5(J),  J=l,  24) 
c 

c... Lower  Quartiles 

read(l,’(24A5)’)  (C5(L).L=1,24) 
write(2,'(/,a,/)’)'  lower  quartiles’ 

write(2,’(12(3x,A5))’)  (C5(J),  J=l,  24) 
c 

c.. .Upper  Deciles 

read(l,’(24A5)’)  (C5(L),L=1,24) 
write(2,'(/,a,/)’)’  upper  deciles’ 

write(2,’(12(3x,A5))’)  (C5(J),  J=l,  24) 
c 

c.. .Lower  Deciles 

readd, '(24A5)’)  (C5(L),L=1,24) 
write(2,’(/,a,/)’)’  lower  deciles’ 

write(2,’(12(3x,A5))’)  (C5(J),  J=l,  24) 
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c 

40  continue 
c 

500  continue 
»  end 
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